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ABSTRACT 

Luminous blue compact galaxies are among the most active galaxies in the local universe in terms of their star formation rate per 
unit mass. They are rare at the current cosmic epoch, but were more abundant in the past and may be seen as the local analogs of 
higher redshift Lyman Break Galaxies. Studies of their kinematics is key to understanding what triggers their unusually active star 
formation. In this work we investigate the kinematics of stars and ionised gas in Haro 11, one of the most luminous blue compact 
galaxies in the local universe. Previous works have indicated that many such galaxies may be triggered by galaxy mergers. We have 
employed Fabry-Perot interferometry, long-slit spectroscopy and Integral Field Unit (IFU) spectroscopy to explore the kinematics 
of Haro 11. We target the near infrared Calcium triplet, and use cross-correlation and penalised pixel fitting techniques to derive the 
stellar velocity field and velocity dispersion. Ionised gas is analysed through emission lines from hydrogen, [Oiii], and [S iii]. When 
spectral resolution and signal to noise allows we investigate the the line profile in detail and identify multiple velocity components 
when present. The spectra reveal a complex velocity field whose components, both stellar and gaseous, we attempt to disentangle. 
We find that to first order, the velocity field and velocity dispersions derived from stars and ionised gas agree. Hence the complexities 
reveal real dynamical disturbances providing further evidence for a merger in Haro 11. Through decomposition of emission lines we 
find evidence for kinematically distinct components, for instance a tidal arm behind the galaxy. The ionised gas velocity field can be 
traced to large galactocentric radii, and shows significant velocity dispersion even far out in the halo. If interpreted as virial motions 
this indicates that Haro 11 may have a mass of ~ lO^^M©. Morphologically and kinematically Haro 11 shows many resemblances with 
the famous Antennae galaxies, but is much denser which is the likely explanation for the higher star formation efficiency in Haro 11. 

Key words, galaxies: evolution - galaxies: kinematics and dynamics - galaxies: individual (Haro 11) - galaxies: starburst - galaxies: 
interactions - galaxies: dwarf 


1. Introduction 

Understanding how galaxies form and evolve is one of the major 
goals of contemporary astrophysics. This involves studying the 
processes which form galaxies’ properties and how they change 
with cosmic time. At high redshift (z) our view of these pro¬ 
cesses is encumbered by poor spatial resolution and low signal- 
to-noise data. Low-z analogues of the objects we see at high red¬ 
shift are therefore, if handled carefully, an indispensable com¬ 
plement to direct studies of the distant universe. 

High-z galaxies selected by the Lyman-brea k technique 
(Lyman-break galaxies, LBGs; [Steidel et al.|1999| ) have proved 
to be important for tracing the star forming galaxy population 
at z > 3. So called Lyman Break Analogues (LBAs) are low-z 
galaxies with gross properties (e.g. UV luminosity and su rface 
brightness) similar to those of LBGs ( [Heckman et al.|20()5] ). The 
number density of LB As is much lower than that of high-z LBGs 
( jHoopes et al.|2007 and references ther ein), and the c losest ex¬ 
amples are VV 114 ([Grimes et^|20Q6| ) and Haro 11 ( [Bergvall[ 
[et al.|2000} [Hayes et al.|2007[ ) both at z ~ 0.02. Haro 11 is one 


of the most luminous Blue Compact Galaxies (BCGs) known 


Send offprint requests to: G. Ostlin 
* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, under observing programmes 71.B-0602, 
074.B-0771(A), 074.B-0802A 


( Kunth & Ostlin 2000{ IB erg vail & Qstlin|[2002 ). Its far UV lu- 
minpsity is Lpuv = 10^^"^L© or 03L*^y '^_^ if compared to the 
z ~ 3 LBG luminosity function ( [Steidel et al.|1999[ ). 

BCGs have attracted much attention due to their properties 
being similar to those expected for young galaxies in the distant 
univ erse: low metallicity, sma ll size, high specific star formation 
rate ( [Searle & Sargent[[1972[ ). BCGs are however an ill defined 
category, including a mixed bag of gross properties ( [Gil de Paz 
[et ar|2Q0^[Micheva et al.|2012a|b[ ). Hence not all BCGs are ex¬ 
treme in terms of metallicity and specific star formation rate, but 
it is among the BCGs that we find the most extreme and efficient 
star f ormation activities in low mass galaxies in the local uni¬ 
verse ( Adamo etal. 2011 ). Whereas the general triggering mech¬ 
anism for starbursts in BCGs is still uncertain, there are indica- 
tions that many of t he luminous BCGs are triggered by mergers 


( Ostlin et al. 2001 Puech et ^[2006 Bergyall & Qstlin[[2002 


[Gumming et al.[[2008[ [Perez-Gallego et al.[[2Qll[ ). One of the 
most pregnant examples is Haro 11, whose morphology is very 
irregular and show many features typical of merging galaxies 
(see Eig.[^. Haro 11 has three bright starburst knots: the south¬ 
western one named ffi, the north-western one !B and the eastern 
one C (see. Eig. and also Vader et aL] (1993); Hayes et al. 
( [2007[ )). In addition, a chain (which we will refer to as the “ear”) 
of young massive star clusters extends between knots and 
Despite its small size and moderate luminosity it forms stars at 


1 













































Ostlin, Marquart, Cumming et al.: Kinematics of Haro 11 


a prodigiou s rate (20 - 30 Mo/yr, |Hayes et"ST]|2007t [Madden 
et al.|2014 ) and with an infrared luminosity of > lO^L© it also 
qualifies as a Infrared Luminous Galaxy (LIRG). 

To further put Haro 11 into context, we compare it to the 
local universe luminosity functions (LFs) of galaxies: Haro 11 

= -18.0 of 

c 


has Lpuv = (compared to M^^y 

et al. 2005), the Ha luminosity is Lna = 


4 4 r* 


Wyder 

▼ 

"Ha 


6.6 X lO'^^erg/s/ cm^ Gallego et al. 1995 1, r-band Lr = 0.9 
(M* = -20.82, [Blanton et al.|2003j ) and = 0.36 L* (M* = 

Tl), whe 


-20.97, [Kochanek et al.[[200lp , where litterature values were 
rescaled to the currently used Hq and photometric system. Since 
the K-band luminosity of Haro 11 is st ill dominated by the b right 
central (r < 5kpc) burst component ( [Micheva et al. 
underlying host galaxy (contributing ~ 30% to the K-\: 


3and 


2010), the 


umi- 


nosity) has a stellar mass of the order of 10% tha t of a typical 
quiescent L* galaxy. The stellar mass estimate of 


Ostlin et al. 


( 200l| and the observed SFR yields a specific SFR of (> 10 




yr or equivalently a birthrate parameter of Z? > 15. 

Studies of LBAs also point to mergers as likely triggers of 
the starbursts that seem to give LBGs their characteris tic proper¬ 
ties ( [Gongalves et al.[20rm Qverzier et al.[2008l[2009| ). In galaxy 
mergers, the detailed dynamical interplay of dark matter, stellar 
components, cold gas inflow, and feedback from triggered star 
formation determine the progress of the starburst and how long it 
can carry on before being quenched, and to what extent its prod¬ 
ucts are made available for chemical enrichment in the merger 
remnant itself or ejected to the intergalactic medium. Such pro- 
cesses have been modelled in large simulations (e.g. [Nagamin^ 
et al.|200^ [Night et al.|2006[ ) but are still poorly resolved, and 
star formation and feedback treated through simplified recipies. 
Local LBG analogues like Haro 11 have the potential to provide 
real observational constraints for such simulations and improve 
our ability to understand galaxy evolution over cosmic time. 

Previous kinematical studies using the Ha line have shown 


that BCGs have very complex kinematics (Ostlin et al. 1999 
20011). Haro 11 presen ts m ultiple kineniatical components (see 


also James et al. 


2013)and 


Ostlin et al. 


( 200lt found that rota¬ 


tion could not support the observed stellar mass. The line width, 
if interpreted as virial motions, indicate a mass of 2 x 10^^ Mq 
consistent with the stellar mass. However, given the very intense 
star formation, it is possible that the ionised gas kinematics in 
BC Gs is driven by feed back and hence do not trace the potential 
(e.g. [Green et al.[20TQ] ). A test of this hypothesis can be provided 
by observations of the kinematics of the stars in BCGs. So far 
results for only a handful B CGs exist in the literature (ESQ 400- 
43, He2-10, ESQ 338-04; Ostlin et al. 2004 Marquart et al.| 
2007t[^mming et al.|2008[ ) but do not show significant differ¬ 
ences be tween the stellar and ionised g as velocity dispersions 
(see also [Kobulnicky & Gebhardt|2000| ). In this paper, we use a 
wide array of spectroscopic data to analyse the kinematical sta¬ 
tus of the stellar and ionised gas components of Haro 11. We list 
some of Haro 11’s basic properties in Table 

The rest of the paper is organised as follows: in Sect. 2 we 
present the observations and reductions from four different ESO 
telescopes and instruments. In Sect. 3 we describe how the stel¬ 
lar kinematics was derived, and in Sect. 4 we do the same for the 
ionised gas kinematics. In sect. 5 we present the results and inter¬ 
pret them, including a comparison with the Antennae galaxies, 
while sect. 6 contains the conclusions. 

Throughout the paper we adopt a luminosity distance of 82 
Mpc (m - M = 34.6) and a scale of 0.4 kpc arcsec"^ valid for 
Hq = 73 km/s/Mpc. 


Table 1. Basic properties of Haro 11 


alias 

ESO350-IG38 

(1) 

coordinates 

00^36”^52^7 -33°33'17" 

(1) 

systemic velocity 

6175 kms"^ 

(1,2) 

distance D 

82 Mpc 

(1) 

metallicity 12+[0/H] 

7.9 

(3) 

Mb 

-20.06 

(3) 

B-V 

0.58 

(3) 

V-K 

2.33 

(3) 

Tpuv 

lO^O-^ Lo 

(7) 

Tfir 

lO^i i Lo 

(7) 

Lua 

3.0 X 10"^^ erg/s/cm^ 

(8) 

Mo-(Ha) 

1.9 X 10^0 Mo 

(4) 


1.6 X 10^0 Mo 

(4) 


_ ^ DOpk (3) Bergvall & Ostlin 

( [2002t , rescaled to Hp = 73, (4)[Ostlin et al. ([2001 1, (7) Ostlin et aTb 


( [20091 ) (8) |Hayes et al.[ ( [2667l > 


2. Observations and reductions 

In this paper we make use of several spectroscopic observa¬ 
tions of Haro 11 all obtained at the ESO telescopes in Chile; 
for a summary see Table In addition we show images from 
the NASA/ESA Hubble Space Telescope (HST), obtained under 
general observer programmes 10575 and 10902 (PI Ostlin). 

2.1. Long-slit spectra with FORS2 

Our long-slit spectra were taken with FORS2 (under ESO pro¬ 
gramme 07 LB-0602) at moderate and low resolution, using two 
slit orientations on the galaxy, see Table and Fig. One, at 
position angle (PA) 57.5° passes through both knot C and knot 
(hereafter referred to as slit C^); the other at PA 140° (here¬ 
after referred to as slit S) passes through kno t tB and the galaxy’ s 
apparent kinematical centre (see Fig. [T] and [Hayes et al.|2007] ). 
We obtained spectra for both position angles using the high res¬ 
olution G1028z grism and a slit width of 0.7", and for slit C7R 
only also using the lower resolution G600B and G600RI grisms. 
The higher-resolution spectra cover ~ 7700-9400A and include 
several transitions of the H i Paschen series, nebular lines such as 
[Sill]9069 and the Can triplet. The lower-resolution data cover 
the spectra region from 3330 A to 8450 A. The relatively narrow 
slit width was chosen to improve spectral resolution {R ~ 4800 
and 2000, for the high and low resolutions grisms respectively). 
A number of template stars (of spectral type F8 to K2), to be 
used in the cross-correlation analysis, were also observed with 
the same setup. The instrument resolution was measured from 
unresolved sky lines. 

The low resolution spectra were obtained and reduced in a 
standard manner with wavelength calibration provided by arc 
lamp exposures. 

The G1028z observations were designed to optimise subtrac¬ 
tion of the strong sky emission from atmospheric OH and O 2 . 
For both slit positions, we took a number of identical exposures 
with the galaxy placed at different positions along the slit. The 
nodding step between two consecutive exposures were set to al¬ 
ways have different length and being larger than the apparent 
size of the galaxy. 

First, a constant background value was subtracted from each 
exposure, with the level determined between the OH lines at 
8510 A. Then adjacent pairs of exposures were subtracted from 
each other, each time scaling the subtracted frame so that resid¬ 
uals from the OH and O 2 lines were minimised. Wavelength cal- 
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Fig. 1. The colour image shows Haro 11 imaged with HST/NICMOS/NIC3 in the F160W filter (~H-band). Based on a 5ks exposure 
obtained under general observer programme 10902. North is up, East to the left. The scale is indicated with a solid white bar. The 
perturbed morphology is evident and demonstrates that the stellar mass distribution of Haro 11 is very asymmetric. Telltale merger 
features like sharp edges and tails are clearly visible. The apertures for the FLAMES/ARGUS and VIMQS/I FU spectra are indi cated 
as grey semi-transparent rectangles. The greyscale image on the upper left shows the HST/ACS Ha image ( Ostlin et al.||2009 ) with 
the three knots !B, C and the chain like structure referred to as the ’ear’ indicated. The positions and orientations of the two 
FORS2 slits are also marked. 


Table 2. Log of spectroscopic observations 


Date 

(UT) 

Seeing 

('') 

Instrument 

Wavelength 
coverage (A) 

Resolution 

(kms"A 

Aperture 

Total integra¬ 
tion time (s) 

1999 09 03 

1.0 

3.6m/CIGALE 

6693-6704 

29 

256x256 pixels (O'.'94 each) 

3180 

2003 08 06 

0.7 

VLT/FORS2 

7675-9420 

63 

slit ^ (through knot !B), PA - 140° 

6075 

2003 08 07 

0.7 

VLT/FORS2 

7675-9420 

63 

slit CJl (through C and Jl), PA = 57.5° 

5400 

2003 08 07 

0.65 

VLT/FORS2 

3330-6120 

160 

slit Cyi, PA = 57.5° 

180 

2003 08 07 

0.65 

VLT/FORS2 

5120-8450 

147 

slit CJA, PA = 57.5° 

180 

2004 11 06 

0.65 

VLT/FLAMES/ARGUS 

8180-9340 

29 

22x14 element IFU (a'52 each), PA 90° 

14400 

2004 10 07 

0.7 

VLT/VIMOSAFU 

4150-6200 

124 

40x40 element IFU (a'66 each), PA 45° 

7100 


ibration was carried out in two dimensions using OH sky lines 
in the non s ubtracted frames, and wavelengths from |Osterbrock 


et al. ( 1996) . To remove pixel-to-pixel variations in the response 


of the CCD, we divided each pair-subtracted frame by a dome 
fiatfield. We found that the OH and O 2 lines varied somewhat 
differently from each other from exposure to exposure. They are. 


however, easy to separate. In the region of interest, O 2 emission 
dominates the background between 8600 and 8720 A, so we 
were able to perform scaling and subtraction for the two cases 
separately, and then spliced the spectra together again after fiat- 
fielding. The spectra were registered and median combined to 
produce the final frames for use in the subsequent analysis. We 
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extracted the spectra three detector rows at a time ((YflS) and data cube as measured by the width of the instrument profile is 
flux-calibrated by comparison with standard stars LTT 7379 and 29 km/s FWHM. 

LTT 7987. 


2.2. Integral field spectra with FLAMES/ARGUS 


Integral held unit (IFU) spectroscopy was taken using the 
FLAMES integral held unit ARGUS, connected to the spec¬ 
trograph GIRAFFE ( Pasquini et al.|2Q02 ) under observing pro¬ 
gramme 074.B-077L We used the 0'.'52 pix“^ scale with grating 
setup LR8 giving a spectral coverage from 8200 to 9380A with 
spectral resolution R ^ 10400. The 22 x 14 lenslets result in a 
fleld-of-view of 1174 x 7'.'3 placed with the long side in east-west 
direction (see Fig.[^. Three spectral template stars (of spectral 
type KOII, KOIII, G8Iab) were observed with the same setup. 

The data were reduced in the same manner as in [Marquart 
et al. (|2007|) using the software and methods of [Piskunov & 
Valenti ( 2002| ) for bias subtraction, flat-fielding, wavelength cali¬ 
bration and optimised extraction of each spectrum. Sky emission 
was subtracted using simultane ously observed spect ra from the 
separate 15 sky fibres (see e.g. [Kaufer et al.||2003| ). As a con¬ 
sistency check, we integrated the spectra over the whole galaxy 
and derived the systemic velocity to be 6179 kms“\ in good 
agreement with previous results (cf. Table [^. 


2.3. Integral field spectra with VIMOS 

Haro 11 was also observed with the VIMOS IFU using the HR- 
blue grism in program 74.B-0802 (PI Bergvall). These observa¬ 
tions, originally designed for a different scientific purpose, were 
obtained at a pointing well outside the main body (15"to the NW, 
see Fig. [^, but contain velocity information on the ionised gas 
that was deemed useful for the current investigation, and there¬ 
fore these data were included in our analysis. 

The VIMOS/IFU with HR-blue grisn|^and scale O'.'67 per 
fibre delivers a held of view of 27" x 27" over the wavelength 
range 4150-6200 A. Of interest for the present work is that we 
detect [Oiii]/I/I4959,5007 and often also Ufi over most of the 
held of view. The spectral resolution is R ^ 2550, verified by 
the width of sky lines in the spectrum. The IFU spectra results 
from a single 7200s integration. In order to remove cosmic rays 
in this single deep exposure, we first normalised each spectrum 
to the continuum near [O iii] and took the median of 5 x 5 pixels, 
constructing a binned spectral image with 3'.'3 resolution. Since 
VIMOS/IFU lack dedicated sky fibers, the sky level was esti¬ 
mated from a dedicated blank sky exposure of 1800s next to the 
target. Some regions of the IFU showed a significantly higher 
residual background level and were not included in the further 
analysis. 


2.4. Fabry-Perot interferometry 

Fabry-Perot interferometry was collected with CIGALE 
( jAmram et al.|199l'] ) mounted on the ESO 3.6-m telescope, un¬ 
der program 63.N-0736, using an etalon with free spectral range 
(FSR) of 390 kms“^ The spatial sampling was 0.94"/pixel. The 
etalon was scanned in 64 channels (each 6.1 kms“^ wide) for 
a total integration time of 3200s. The data was reduced with the 
Adhoc ( Boulesteix|1993 ) software, and a Gaussian smoothing of 
T.'5 was applied. The effective spectral resolution in the reduced 


^ Note that this is the old HR-blue grism. The new HR-blue grism 
has lower spectral resolution and a different wavelength range. 


3. Stellar kinematics from the Caicium tripiet 

Observing the stellar velocity held in metal-poor starbursting 
galaxies is very challenging due to intrinsically weak stellar ab¬ 
sorption lines and contamination from ionised gas emission. In 
galaxies like Haro 11 with very strong emission line spectrum, 
the Baimer lines become completely emission dominated. The 
Calcium triplet, which originates from photospheric absorption 
in red giants and red super giant stars, is one of very few prac¬ 
tically accessible probes. However the Calcium triplet lines are 
partly blended with Paschen lines which needs to be dealt with 
before stellar velocities can be derived. 


3.1. Removing Paschen emission lines 


We subtracted the lines of the Hi Paschen series from our 
FORS2 spectra according to the scheme presented in jCummm^ 
et al. ( 2008 ). As a template for the higher order Paschen lines, 
we constructed a model profile by adding the profiles of the three 
strongest H i lines in our spectra that are not blended with the 
Calcium triplet lines (Pa 10, 11 and 12). Where the signal to 
noise ratio was low, we used [S iii] 49069 as a template. 

We adjusted the width, flux and velocity of the model lines 
in order to give the best subtraction around those Paschen lines 
which do not blend with Can absorption. Finally we masked the 
few other emission lines seen in the spectra, O i 48446, [Cl ii] 
48579 and Fe ii 48617. 

For the ARGUS spectra, we followed the same procedure to 
subtract Paschen emission lines, using Pa 10 as the sole profile 
template due to the limited wavelength coverage. 


3.2. Stellar kinematics from FORS2 slit spectra 

We derived the stellar kinematics by cross-correlation of the 
cleaned spectra using FXCOR (based onjTonry & Davis||1979 ) 
and the same set of templat e stars as in j Cumnnng et al. |( |2008| ). 
We followed the method of Ostlin et d.| |2007| ) for deriving the 
resulting errors in the velocity dispersion. Figure shows how 
the cross correlation function (CCF) peak and the [Siii] pro¬ 
files vary along each slit. In the right panels of figures and 
we show the resulting stellar velocities (relative to the sys¬ 
temic velocity 6175 kms“^) and velocity dispersions (with or¬ 
ange squares plus error bars). Independently, we derived the ve¬ 
locities and dispersions using the penalised pixel-fitting (pPXF) 
method of jCappellari & Emsellemj ( |2004| ) (shown as magenta 
circles in Fig. an W 

The derived velocities and velocity dispersions are subject 
to systematic errors from the Paschen line subtraction, in 
particular where the Paschen emission is strong and the Calcium 
triplet weak. In Fig those CCFs for which the peak amplitude 
was below 0.25, or the Paschen-subtraction was not deemed 
reliable, are shown in grey, and these points were discarded 
from the subsequent analysis. To quantify the effects of over- 
and under-subtraction of the Paschen lines, we additionally 
investigated two extreme cases: 


(i) We increased the Paschen-line intensities as much as possi¬ 
ble without inducing artificial absorption features , all the while 
adjusting the velocity dispersion and velocity of the nebular line 
model to allow the maximum intensity to be subtracted. 
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5500 6000 6500 7000 


5500 6000 6500 7000 




Fig. 2. Variation along the slit of the cross-correlation function (CCF, first and third columns) from fxcor and the [S iii] A9Q69 line 
profile (second and fourth columns). The leftmost two columns show the results for slit B and the rightmost columns the results foj 
slit Where the CCF peak value is smaller than 0.25 or the Paschen-subtraction is deemed unreliable, the CCF is shown in grey. 
The emission-line profiles have been continuum-subtracted. Each sub-panel is labelled with the position of the extraction window, 
where the origin is set to the position where the two slits cross (see Fig.[^. Knot B is at position -rT.'9 along slit and knots C and 
are at -2'.'2 and -f 3" along slit C^, respectively. 
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(ii) In the other case we instead kept the subtracted intensity as 
low as possible without leaving any residual emission features at 
the unencumbered Paschen lines, again letting the velocity and 
velocity dispersion take any values to assure this. 

The spectra resulting from the two alternative subtractions 
were then cross-correlated as before. The results bracket the pos¬ 
sible range of velocities allowed by the data, and give a good 
estimate of the maximum systematic errors. In figures and 
we show, with upwards and downwards triangles, the velocities 
and dispersions resulting from these extreme Paschen subtrac¬ 
tions. The comparison with the pPXF results give further infor¬ 
mation on the systematics. Sometimes the differences are larger 
than those represented by the error bars, and other times smaller. 
In the latter case the obtained result should be regarded as quite 
robust. The stellar velocities are found to be insensitive to the 
template star used in the cross-correlation, so we treat the data 
for the bright KO giant HD 164349 to be representative. 

For slit we see that the stellar velocities and dispersions 
resulting from the two different methods (FXCOR and PPXF) 
and extreme Paschen-subtractions agree very well, and hence 
the overall velocity and dispersion profile must be regarded as 
very reliable. For slit ^ the agreement is worse, partly due to the 
emission lines being stronger along this slit. 


3.3. Stellar kinematics from FLAMES/ARGUS 


A fter Paschen-subt r action we followed the procedure described 
in [Marquart et al.| ( |2QQ7| ). However, poor signal-to-noise and 
strong Pa emission lines in many of the spaxels meant that we 
could not derive a detailed spatial map of the stellar velocities. 
After binning 2x2 we could determine the stellar velocity in the 
region around C where the Paschen emission is relatively weak, 
finding good agreement with the FORS2 result. The only region 
where we get reliable velocity data and that was not covered by 
by FORS2 is 2" South-East of knot C. Here the stellar velocities 
are the same as those for the gas (derived from [S iii]). We will 
not discuss the stellar kinematics from ARGUS further. 


4. Ionised gas kinematics 


Velocities for the ionised gas are calculated using rest wave¬ 
lengths from [van Hoof| (|1999|) and (for [S iii] only) |Hippelein| 
|& Muench|([T98T]): 


4.1. Gas kinematics along the FORS2 slits 

The analysis of the FORS2 ionized gas velocities and their er¬ 
rors follow ( [Cumming et al. 2008[ ). In general we use [Siii] 
49069, the strongest line in the G1028z spectra. Along slit po¬ 
sition CjTI we also observed Ra with grism G600RI, which al¬ 
low us to probe the velocities to somewhat larger radii than for 
[Sill], albeit with lower spectral resolution. The [Siii] emission 
line profiles are shown in Fig and it is obvious that they are 
non-Gaussian at many places. Since the ARGUS data has higher 
spectral resolution we restrict our analysis of the detailed line 
shape to this data, see next subsection. In Fig. and we show 
the velocities and velocity dispersions for stars and gas derived 
from the assumption of a single Gaussian component. We note 
that the regions coinciding with knots JTI, S and C have smaller 
velocity dispersion than the average of their surroundings. 


4.2. ARGUS 

We used [S iii] 49069, the strongest emission line in our spec¬ 
tra from ARGUS, to derive spatially resolved kinematics of the 
ionised interstellar medium. The top row of Fig. [^displays maps 
of amplitude, velocity and velocity dispersion resulting from a 
fit to the line using a single velocity component. Because the 
line shape at this spectral resolution is non-Gaussian over large 
regions of the field-of-view, we allowed for two additional de¬ 
grees of freedom in the form of the parameters and /z 4 from the 
Gauss-Hermite polynomials. They represent skewness and kur- 
tosis (spikiness/boxiness) respectively and generally allow the 
fitted velocity to better correspond to the peak of the line profile. 

The [S III] line shows a complex kinematic structure and sev¬ 
eral individual line profiles indicate the presence of multiple ve¬ 
locity components. We next decomposed the [S iii] line profiles 
into several components by fitting multiple Gaussians in the re¬ 
gions where it is clear that the profile cannot be described by 
a single component. We use the smallest numbers of compo¬ 
nents deemed necessary from the shape of the line profile and 
the quality of the data. This means that for some pixels we have 
one component, for some two and for others three. Even three 
components is sometimes insufficient (e.g. lower right panel la¬ 
belled ’4’ in Fig.[^ but we refrain from fitting more that 3 profile 
due to the ambiguity in the interpretation. 

The line decomposition in some places reveal up to three dy¬ 
namical components - how can we spatially make sense of these 
components? There is a principle problem of multi-component 
fits in that one needs to arrive at a spatially consistent identifi¬ 
cation and labelling. The question whether a component at one 
place in the galaxy is related to a particular component in a 
nearby region is non-trivial as each component may vary spa¬ 
tially both in velocity, width and amplitude. We use spatial con¬ 
tinuity over several pixels as a criterion for manually deciding on 
how to sort the components in an iterative fashion. This means 
that we tried to avoid discrete jumps in velocity and line width 
for a given component. Rows 2-4 in Fig. [^display maps of the 
resulting kinematical components: the first component is wide 
spread over the field. The second component forms a kinemati¬ 
cally cold structure that is well defined in space, intensity, veloc¬ 
ity and line width and which we associate with the ’’ear”. The 
third component is somewhat more elusive. The first compo¬ 
nent represents everything that does not belong to the second 
or third component or when the line shape did not obviously re¬ 
quire more than one component. 

In addition to the sorting problem, there is another uncer- 
taintiy that affect the velocities and dispersions shown in Fig.|^ 
Multiple components can only be clearly identified if they have 
different amplitude or are separated in velocity by an amount 
comparable to the FWHM of the individual lines. This means 
that even where a single component may seem sufficient, there 
might in reality be several distinct regions projected along the 
line-of-sight. If the velocity separation and relative intensity of 
different components change, at some positions they may blend 
and two relatively narrow slightly separated components may be 
mistaken for a single broad one, leading to velocity and cr jumps. 


4.3. CIGALE 

The Hof velocity field from CIGALE was analysed by fitting 
a Gauss-Hermite polynomial to the line profiles in each pixel. 
The resulting velocity field (hi) and velocity dispersion (/z 2 ) are 
shown in Figure]^ There is a good agreement with the [S iii] ve¬ 
locity field (Fig.^. In order to allow comparison with the results 
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emission lines 


Ca II (template HD 164349, KO II) 



Fig. 3. Kinematical results along FORS2 slit B showing (from top to bottom) the flux level, the line of sight velocity dispersion, and 
the heliocentric radial velocity (relative to the systemic velocity of 6175 km s“^). Zero on the abscissa is deflned as the point where 
our two slits cross. The left panels show results for the [Siii] A9069 (red open circles), and Hi Paschen emission lines (weighted 
means of all measured lines; blue fllled circles). The flux scale show the normalised strengths of the emission lines. The dotted 
vertical lines at 2" and 5" marks the position along the slit of knot !B and where the slit crosses the loop of structures referred to 
as the ’ear’ (see also Figs. and [^, respectively. The dotted horisontal line in the cr-plot is arbitrarily put at 75 kms“^ The right 
panel shows at the top the spatial variation of the normalised continuum flux near the Calcium triplet (orange), in arbitrary units. 
The red solid line repeats the [S iii] flux distribution from the left panel. The middle and bottom panels show the derived stellar 
velocity dispersion and radial velocity, respectively. Again, the solid red line repeats the [S iii] results from the left panels. The 
results from cross-correlation is shown as orange squ ares and the pPXF results as magenta open circles. The grey triangles represent 
the extreme values derived (as described in sect. |3.2| ) from the cross-correlation method when systematics in the Paschen emission 
line subtraction is taken into account. At the distance of Haro 11, 10" corresponds to 4 kpc. 


of 


Ostlin et al. 


(1999 20011 we have also extracted a rotation 


curve alo ng the major axis defl ned by the new photometry pre¬ 
sented in |Micheva et al.| ( 2010| ), where the shape of the K = 23 
mag/arcsec^ isophote was fitted with a position angle PA = 103° 
and an inclination i = 47° (see Fig.[^. The line of sight Ha line 
widths are cr ^ 60 kms“^ at and C, and ~ 80 kms“^ at S. 
The broadest lines are seen in the ’overlap’ region (i.e. approxi¬ 
mately the point where the two FORS2 slits cross and where the 
ARGUS data reveals at least 3 velocity components in the [S iii]- 
line), and also 4" south of this region and in addition 2-3" north 
and east of C. 


4.4. VIMOS 

In the VIMOS data, the H/5, [Oiii] 44959, and 45007 lines are 
detected in all binned pixels. We simultaneously fitted the [O iii] 


lines with a single redshifted Gaussian each, connected by the 
intensity ratio 1:3, a common width and a separation of A4 = 
(5007 - 4959)(1 -f v^/c) A where is the fitted radial velocity. 
We independently fitted the HJS line with a single gaussian. 

The fleld-of-view of the VIMOS IFU was placed at a 45° an¬ 
gle north-west of knot S. Since we are interested in how the gas 
velocity and line width changes with distance from the starburst, 
we collapsed the velocity held along the direction of slit B, re¬ 
sulting in a single set of velocity versus radius data, shown as 
black dots in Fig. [7] Scattered light could potentially be an issue 
with VIMOS IFU data ( |Monreal-Ibero et al.|2006| ), however, our 
pointing is far from the bright centre which is well outside the 
held of view of VIMOS/IFU. 

The HJS data is quite noisy but on average the [O iii]/HyS ra¬ 
tio is ~ 5, similar to what is seen in the centre (Bergvall & 
Ostlin 2002). We And that the velocities to match those of 
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Fig. 4. Same as Figurebut for slit CJK, and with the addition of Ha measurements (cyan triangles). In the right panels, the cyan 
full drawn line repeats the Ha measurements from the left panels, to allow for comparison of the stellar and ionised gas kinematics. 
Dotted vertical lines mark the positions of knots C (left dotted line) and (two rightmost dotted lines), see Fig.[^ 


the [O III] lines within the uncertainties. Therefore we only fur¬ 
ther investigate the latter. The [Oiii] lines show no apparent 
structure, such as asymmetries or double components. They are 
however, significantly broader than the instrument profile, which 
we measured to be = 124km s“^ (or equivalently 

t^instr = 53km s“^) . In the high signal to noise bins, the raw 
width as determined from a Gaussian is CTraw = 100+10 km s“\ 
implying an instrument corrected line width of cr 85 + 10 
km s~^ i.e. very similar to the average value in the centre of 81 
km s“^( Ostlin et al. |2001| ), see also Fig. There is a weak de¬ 
cline in (T with increasing radius. 


5. Results and analysis 

Figures and show how the velocity, a and fiux vary along 
the FORS2 slits, for selected emission lines (left panels) and for 
Can (right panels). The results for slit CJi (Fig.[^ show a clear 
correspondence between the velocities and dispersions for gas 
and stars. The velocity curves are not identical, but overall, stars 
and gas have very similar kinematical properties. Slit B adheres 
to this trend even if the poorer data makes it less obvious. 

The ARGUS ionised gas velocity field on closer inspection 
showed evidence for multiple kinematical components, as dis¬ 
cussed below. Also the stellar velocities may be multi compo¬ 
nent in nature as hinted by the shape of some of the CCFs in Fig. 


1^ but the current data does not allow us to reliably investigate 
this further. 

The similarity of the stellar and gaseous velocities sug¬ 
gest that have a similar physical origin, namely virial motions. 
Whereas gas is known to be affected by feedback, this appar¬ 
ently has not affected the line core, perhaps because gas that is 
affected by feedback tends to become diluted and have low in¬ 
tensity. The motions of stars are, on the other hand, not affected 
by feedback. 

5.1. Analysis of the [Sin] velocity field 

The high spectral resolution of the ARGUS IFU-spectra allows 
us to have a closer look at the shape of a gaseous emission line. 
If one measures a line shape that is clearly the sum of multi¬ 
ple components with different velocities, it is safe to assume that 
they originate from physically distinct regions since it is hard 
to imagine a single region with two or more distinct bulk mo¬ 
tions. This means that we have a handle on the third dimension 
of Haro 11, perpendicular to the plane of the sky. 

Fig 1^ shows the result and some examples of fitted line pro¬ 
files. Several results emerge from this: 

(1) There is a main velocity component in the sense that it 
can be seen across the galaxy. There seems to be no indication of 
a discontinuity as one would expect in a system that has not yet 
merged. The decomposed line has its broadest (cr ^ 90km s“^) 
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Fig. 5. Maps of the field-of-view of FLAMES/ARGUS showing the derived properties of the [S iii] A9069 line. The tick marks have 
a spacing of 1". The top row shows the results from a single component fit, from left to right: intensity in arbitrary logarithmic 
units (with the location of the 3 knots indicated with crosses), velocity and velocity dispersion. Since we find that over large areas, 
a single velocity component is insufficient to reproduce the line profile, we present an attempt to decompose the line by fitting 
multiple gaussians. The number of components was chosen conservatively and additional components were only included in when 
deemed necessary to produce a decent fit. The results are presented in the three middle rows. The bottom row shows examples of 
measured [S iii] line profiles (black solid line), and our multi-component fits to them (first: blue dotted, second: green dash-dotted, 
third: yellow dashed). The red solid line is the sum of the components. The x-axis displays velocity, relative to the systemic velocity, 
in km s“^. The positions that correspond to the spectra are shown by numbers over-plotted on the panels above them. 
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Fig. 6. Results from CIGALE Fabry-Perot Ha observations of Haro 11. Upper left: Line-of-sight velocity (relative to the systemic 
velocity) with the scale (in kms“^) given by the colourbar beneath. The spatial scale is indicated in the corner and a single Ha 
flux contour (at ~ 5 • 10“^^erg s“^cm“^arcsec“^) is over-plotted for reference. The solid black line shows the line along which the 
rotation curve (in the lowe left panel) was extracted and the dashed lines indicate the +40 deg wedge of points included. Upper 
right: Velocity dispersion map. The same contour as in the first panel is over-plotted. The white circles mark the positions of knots 
tBmd C, and diamond marks the “overlap region”. Lower left: Rotation curve, derived by deprojecting the apparent position of 
each velocity point to its radius in the inclined plane, using the position angle and wedge indicated in the upper left panel. The error 
bars correspond to the spread of the velocity measurements within each one arcsecond wide radial bin. Lower right: Example of Ha 
profiles from the locations of the three knots and the “overlap region”. 


component ~ 1" south-west of knot S (see spectrum ’2’ in 
Fig|?|. This is also the apparently most dusty region (see Fig. 
1^. It is possible that the broad component giving rise to the high 
line width is in fact due to unresolved narrower components. 

(2) The second velocity component show a near spatial co¬ 
incidence with the ’ear’ (Fig. [^. It forms a high velocity com¬ 
ponent that is narrow in space and velocity (cr ~ 30 kms“^). It 
is likely a likely tidal arm, that was flung out in the merger, and 
which may be associated with the ’ear’. 


In a less detailed investigation, the high velocities at the lo¬ 
cation of the “ear” could have been misinterpreted as due to rota¬ 
tion. Structures like this one could easily be hiding in low spec¬ 
tral resolution studies of BCGs or merging galaxies, and high¬ 
light a potential pitfall in interpreting the kinematics of complex 
galaxies through low resolution spectra. 

The kinematical information from VIMOS (see Fig. [ 7 ]) sup¬ 
port this interpretation. The velocities measured further to the 
North-West (Fig.[^ do not coincide with the high velocity of the 
second component but instead joins smoothly with the first com- 
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Fig. 7. Line-of-sight velocity of gaseous emission lines along 
PA -45°. Black dots are the [Oiii] velocities and line of sight 
cr from our VIMOS data. The errorbars for all points represent 
the spread (Icr) of the individual measurements within the radial 
bin, not the measurement error. Green, red and blue dots cor¬ 
respond to the first, second and third velocity component from 
the ARGUS analysis of the [S iii] line (cf. Fig. [^, as derived 
by a cut through the velocity field in the north-west direction. 
Small symbols represent individual spaniels and large symbols 
radial averages with error bars show the spread. The scale is 0.4 
kpc/arcsec. 



solved if not viewed with the HST or adaptive optics. Using the 
baolab package ( |Larsen|1999| ) and the HST/ACS/HRC/F814W 
image, we find the light distribution to be best fitted with a King 
profile with effective radius of reff = 40 + 10 pc. This indicates a 
virial mass for knor C of: 


Mv 


2.4 X 10' 


f ^ 

\50kms V \ 


^eff 

40 pc 


Mo. 


( 1 ) 


Adamo et al. ( 2010| ) used SED fitting to estimate an order 


of magnitude lower stellar mass which however was based on 
the assumption that C was a point source and therefore grossly 
underestimated the total flux. 

In HST images, knot ^ is resolved into a group of bright star 
clusters. The size of the group is reff ^ 200pc and the velocity 
dispersion cr 50 kms“^. Hence the dynamical mass of the 
complex making up knot ^ is Alyir ~ lO^M©. 

For knot tB we measure cr 80km s“^ although it is more 
uncertain than for and C due to the very bright emission lines. 
With a size of rgff ^ 140 pc the mass is of the order ~ 2 x IO^Mq. 

The velocity dispersion (both for gas and stars) have local 
minima at the positions of knots and bright clusters. This is un¬ 
derstandable since in such cases the light is dominated by the 
knot in question and the measured line width is strongly biased 
by the local conditions in the knot, whereas in a random position, 
the line profile has (often unresolved) contributions from many 
velocity components along the line of sight. It does not imply 
that the knots are necessarily dynamically cooler. 

Comparing the stellar and nebular gas velocity field it is 
striking how similar they are on average. While we do not have 
the same possibility to resolve the stellar velocity field into mul¬ 
tiple components, the similarities indicate that Can and nebular 
lines trace the dynamics of the system in a similar way. Stars are 
not subject to outflows and hence, while there may well be out¬ 
flow components present in the nebular spectrum at many points, 
they do not dominate the local fluxes. 


ponent. Hence, the “ear” is a spatially and kinematically distinct 
region. Figure shows that the Ha emission in t he ear is dom¬ 
inated by a chain of compact young star clusters, [Adamo et al.| 
( |2010|). Similar structures as the ear are seen in the Antennae 
(see Whitmore et al.|2QfO| and also Sect. |5.5| in this paper). 

(3) There is a third lower velocity component, mostly around 
knot A and the ’’overlap” region (position 3 in Fig|^. 

The Ha velocities probed by Cigale (Fig.[^ in general agree 
well with those from [Siii] although the spatial resolution is 
worse. The Ha line is markedly non Gaussian over most of the 
galaxy, and is broadest in the overlap region. 


5.2. Stellar kinematics and virial masses of knots Ji,!B and C 

The stellar velocities along slit CJi shows a gradient of about 50 
kms“^ kpc“^ (Fig.|^. The stellar velocity dispersion ranges be¬ 
tween 60 and 150 kms“^ and largely follows that of the gas, 
peaking halfway between knots ^ and C (close to 0"on the 
abscissa, i.e. the overlap region). At this point the ionised gas 
spectrum from all instruments (CIGALE, ARGUS, FORS2) is 
markedly multi-component. Also the CCF structure is sugges¬ 
tive of multiple stellar components, and this is the most likely 
explanation for the very high stellar velocity dispersion found. 

Both the stellar continuum and the absorption line equivalent 
width of the Ca ii triplet are strongest at the position of knot C 
where the velocity dispersion of gas and stars show a local min¬ 
imum of (T 50 + 10 km s“^. Knot C is very compact and unre- 


5.3. The VIMOS far out velocity field 


How can we understand the VIMOS results showing a fiat ve¬ 
locity curve and an emission line width of aios ~ 80 km s“^ at 
galactocentric distances >10 kpc? First we note that this value is 
similar to the average found for a one component fit to the neb¬ 
ular gas (and the stars) in the more central regions. The spectral 
resolution of the VIMOS spectra and the S/N do not allow for 
a multi component analysis. Below we discuss a few alternative 
explanations for the broad lines. 

Thermal broadening: Thermal broadening can be ruled out 
since the implied temperature (~ lO^K) would be unrealistic for 
a gas showing [O iii] in emission. The [O iii] electron tempera- 
ture in the centre (where the line is equally broad) is ~ lO'^K 
( Bergvall & Qstlin||2002 ). Hence, the line broadening must be 
macroscopic and due to bulk motions. 

Rotational broadening: If the far reaches of Haro 11 were 
a rotating disk, the line profile would be broadened by the fact 
that the rotational velocity will have different sign and amplitude 
along the line of sight. Assume we view the galaxy at some im¬ 
pact parameter h according to Fig. and that the galaxy has a 
constant rotational speed Vrot which we observe at an inclination 
i. The observed velocity is Vobs = Vrot sin(/) and from Fig. we 
have that Vobs = 50 kms“^. At Z? » 0 the line would be made 
up by components with velocities ranging from ~ Vobs! V2 ~ 
35km s“^ (at I = Imin and Imax) to Vobs at / = 0, and the induced 
line broadening would be CTobs ~ (Vobs - Vobs/ V2)/2.35 6 
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locity dispersion. Moreover, the mechanical energy output of 
an instantaneous starburst peaks at about 20 Myr (Starburst99, 
Leitherer et al.||1999|). From analy sis of the luminous star clus- 


Fig . 8. Schematic illustration of relevance for the discussion in 
|5.3[ An observer sees the galaxy at an impact parameter h along 
the line of sight (dashed line). The line of sight intersects the 
hypothetical galaxy between Imin and Imax 


km s“^ (where the factor 1/2.35 comes from the transformation 
of FWHM to (t). Hence, rotation give a negligible contribution 
to the line width. 

A radial outflow: A scenario which is compatible with a flat 
velocity gradient and broad emission lines is a uniform radial ex¬ 
pansion, since this does not affect the velocity centroid (i.e. the 
rotation curve) but broadens the line. A spectrum at a projected 
galactocentric distance h would sample velocities with line of 

sight components vios = Vexp X // from to l^ax, 

see Fig At / = 0 the expansion is purely transversal, whereas 
I < 0 and I > 0 produce blue- and red-shifted components, re¬ 
spectively, which broaden the line. The resulting broadening in 


FWHM would be approximately Vexp which hence would have 
to be ~ 200 km s“^ in order to explain the width of the line. The 
line width would be expected to get narrower at larger radii and 
there is indeed a weak trend in that direction. 

However, a real outflow is unlikely to be spherically sym¬ 
metric. For instance pohnson etaL] ( 2012| ) investigated the nearby 
dwarf galaxy NGC 1569 in ionised gas and Hi finding very ex¬ 
tended and irregular bipolar velocity structures. If the far out 
velocity held in Haro 11 is due to outflows, it would be coinci¬ 
dence th^he VIMOS FOV happened to capture them. |Kunth 
|et al.| ( p998 ) found evidence for an outflow (of neutral gas) with 
58 km/s from HST/GHRS spectroscopy in a pointing towards 
the centre of Haro 11. Additional si ght lines were probed by Na 
absorption in |Sandberg et al. ( 2013| ) who found an outflow of 44 
km s“^ towards B but an inflow of 27 km s“^ towards C. In con¬ 
clusion, there is no evidence for a galaxy wide radial outflow of 
the magnitude needed to explain the halo line width. 

Turbulence from feedback: It has been debated whether 
emission line widths are due to virial motions or due to turbu- 
lence created through feedback fr om star formation ([Terlevich & 


ter population I Adamo et al.| ( |2010| ) found that the star formation 
rate in Haro 11 started to increase rapidly 20 Myr ago. Hence, 
the current starburst has had limited time to transport mechan¬ 
ical energy out into the halo. For a scale of 10 kpc and a time 
scale of 10 Myr, the implied velocity is 1000 km/s. We see no 
indications of a wind at this speed, and hence the line width is 
unlikely to have been produced purely by feedback. 

Virial motions: We are left with the most plausible expla¬ 
nation being that the line width is mainly due to macroscopic 
bulk motions of the ions. We can then use the clouds as trac¬ 
ers of the gravitational potential and equation (1) to get an or¬ 
der of magnitude dynamical mass estimate of ~ lO^^M©. The 
rotational velocity (~ 50/sin(/) kms“^at r = 12kpc) provides 
another 5 x 10^/ sin^(/) Mobut is negligible in this context, and 
the ’halo’ of Haro 11 would be dominated by velocity disper¬ 
sion. This assumes dynamical equilibrium and isotropic velocity 
dispersion, both which are questionable for a merging system 
like Haro 11, but the order of magnitude dynamical mass should 
nevertheless be lO^^Moin this scenario. If the merging galaxies 
came in on parabolic orbits, pre virialisation velocities will be 
^^2 higher and the inferred virial mass a factor of 2 too high. 

To conclude, Haro 11 is rotating to the largest scales we can 
probe, and the rotationally supported mass is, under the assump¬ 
tion of circular motions, ~ 10^^ M©. The emission lines in the 
halo region are broad, cr ~ 80 km s“^ We And that the line width 
is likely due to bulk motions indicating a mass of ~ 10^^ M©. 


5.4. Stellar and dynamical masses 

From multi-wavelength data it is possible to estimate the mass 
in stars through comparison with spectral evolutionary synthesis 
models, under the assumption of a stellar initial mass function 
(IMF) and a certain star formation history. This is commonly re¬ 
ferred to as stellar mass. While the resulting age and star forma¬ 
tion history is subject to degeneracies, th e stellar mass is o ften 
mo re robust (for reason s discussed in e.g. Ostlin et al.|2Q01 1. 
Ostlin et ^ ( 2001 ) found that for radii r > 1 kpc, the stellar 


mass distribution in Haro 11 could not be supported by rotation. 
Only by assuming that the observed Ha velocity dispersion trace 
mass could the dynamical and stellar masses be reconciled. Then 
the inferred dynamical mass is 1.9 x 10^^ M© for an effective 
radius of 2.8 kpc and cr = 81 km s“^. 

In this paper we also make use of the stellar popu lation mod¬ 
elling of Haro 11 performed in [Hayes et al. ( |2007| ). In that pa¬ 
per, the stellar population was modelled in a spatially resolved 
fashion using HST multi-band imaging from 1500 A to the I- 
band, and with 2 populations (one young and one old) plus one 
ionised gas component. The underlying assumptions were that 
both populations are single stellar populations, calculated with 
Starburst99 and Geneva tracks with metallicity Z = 0.004 and a 
Salpeter| ( [l955| ) IMF. If we ntegrate the stellar mass from the ref- 


Melnickj 1981 [[Green et al.|20~T0}[Moiseev & Lozinskaya|2012| ). 

The fact that stellar and gas velocity dispersions tend to agree in 
Haro 11 argues against turbulent feedback driving the gas ve- 


erence position where the slits cross (r = 0) and out to r = 10'" 
we And a total stellar mass of 2.4x 10^^ M©. Half the quoted mass 
is contributed by the inner r < 3" that encompasses the three 
main knots. The stellar mass estimate is primarily sensitive to 
the assumed low mass slope and cutoff of the IMF. For more re- 
alistic IMFs wh ere the slope flattens a low mass (e.g. [Scalo 1986[ 
[Chabrier|2Q03[ ) this estimate should be reduced to about 60% of 
the quoted value, or 1.4 x lO^^Mn. The se estimates are of the 
same order as those by Ostlin et al. ( 2001 ) (M* = 1.6x10^^ M©) 
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and Madden et al. ( |2Q14 ) (M* = 1.7 x 10^^ M©), both studies 
using independent data. 

We see a sharp velocity gradient along slit from -6" to 
+4" (i.e. a total extent of 4 kpc) over which the velocity varies 
monotonically with 200 kms“^ Interpreting this in terms of ro¬ 
tation gives a rotationally supported mass of 4.6 X 10^/ sin^(0 
Mq, which for i > 30° is smaller than the mass inferred from the 
line width (if using the same radius). The inclination is not well 
constrained, but rotation likely gives a significant contribution 
in the central (r < 10") region. Further out, we can follow the 
Ha kinematics along slit CJi from -12" to +16" (a total scale 
of 11 kpc). The velocities changes slope on both sides, but the 
gradients are about equally steep, and the same is observed in 
the CIGALE velocity field (Fig.[^. 

Hence the velocity field in this ongoing merger is complex. 
If the velocity gradients are interpreted in terms of gravitational 
motions, there seems to be ample dynamical mass also outside 
the central 4 kpc. The velocity amplitudes and cr are the same 
as in the centre but the spatial extent almost 3 times larger. The 
total dynamical mass for the Haro 11 system may under this as¬ 
sumption approach lO^^M©. This is close to the value inferred 
from the VIMOS far out velocity field discussed in the previous 
subsection 

5.5. Haro 11 - a denser version of the Antennae 

Haro 11 presents a very rich and complex velocity struc¬ 
ture. Given the morphological features seen in the 1.6yum 
HST/NICMOS image (Fig. this comes as no surprise. 
However, Haro 11 presents yet more intriguing features. 
Comparing it to the Antennae galaxies reveal striking similar¬ 
ities. We demonstrate this in Fig. where the HST images of 
the Antenna^ have been rotated and scaled to the same appar¬ 
ent orientation and size as Haro 11. On detailed inspection they 
almost almost look like copies and bright irregular features and 
dust lanes in Haro 11 have close counterparts at the same appar¬ 
ent place in the Antennae and vice versa. 

In Fig. 1^ we highlight, using roman numerals, some con¬ 
spicuous features which are also discussed below: Both galaxies 
present distinct edges in the underlying stellar distribution (I). 
Dust streaks that in the reference system of Haro 11 runs hori¬ 
zontally north of knot C (II). A dusty Hcr-bright region at !B (HI). 
The region which in Haro 11 is called is bright in Ha and has 
many optically bright and blue star clusters SE of it (IV). A faint 
tidal feature south-east of (V). An optically bright compact 
nucleus C with dust lanes leading into it (VI). Also the ’ear’ 
in Haro 11 may have its analogue in the Antennae system (see 
e.g. Sect. 8 in Whitmore et al.|20T0| ). Of course, there are dif¬ 
ferences too: knot ^ is a complex of star clusters rather than a 
nucleus as in the Antennae, and Haro 11 lack extended tidal tails 
(see below). Nevertheless, the striking similarity of the the two 
systems suggest a similar origin - the mergers appear equally ad¬ 
vanced and has presumably been created by progenitor galaxies 
that are similar to both systems and that has had similar orbits 
with respect to each other. However, despite the striking simi¬ 
larities there is one important difference: Haro 11 is physically 4 
times smaller than the Antennae. 

The ionised gas velocity field of the Antennae was studied 
by lAnn’am et al. ( |1992| ) where they used CIGALE, i.e. the same 
Fabry-Perot interferometer as in this work, to study the Ha line. 
The smaller distance and larger physical size allow for greater 


detail to be seen in the Antennae, so we compare gross proper¬ 
ties (in the reference system of Haro 11). The region around C 
in both galaxies show low velocity, the region around have 
significantly higher velocities, and with !B being intermediate. 
Hence, not only the image projected on the sky plane is simi- 
lar in b oth galaxies, but also the relative motions. Amram et al. 
( |1992| ) showed a position-velocity (PV) diagram, extracted along 
the line joining the two nuclei, i.e. corresponding slit C^. We ex¬ 
tracted a PV diagram from our CIGALE data cube along this di¬ 
rection and compare it with that of Antennae in Fig.[^ Again, the 
nearer distance to Antennae allows for more detail and a larger 
number of velocity points to be extracted but overall the two 
PV diagrams are strikingly similar: they share the same overall 
shape and have the same velocity amplitudes (Av 250 km s“^). 
The recent evolution of the SFR in the two galaxies also appears 
to be si milar. From studying the ages and masses of young star 
clusters |Adamo et al.| ( |2010| ) find that the SFR in Haro 11 has 
increased with a factor of 50 over the last 40 Myr. This corre¬ 
sponds to the estimated time of the 2nd pericenter passage of 
the merging galaxies in the Antennae after which the SFR has 
increa sed with a factor of 40 ( [Karl et al.|2010t [Whitmore et al.| 
2007 1 ). Hence the time scales for the two mergers seems to be 
similar. 

While knots and C are separated by 2 kpc, and the nuclei 
in Antennae by 8 kpc (indicating a factor of 4 difference in size) 
the PV diagram in Haro 11 appears relatively more stretched. 
The distance between the extreme velocity values is 3.7 kpc in 
Haro Hand 7.6 in Antennae. This is because in Haro 11 C and 
does not lie at the extreme values, but the velocity curve con¬ 
tinues with the same shape beyond them. This is obvious also 
in Figures and 1^ Hence Haro 11 is more compact in position 
space than in velocity space. Taking into account that the beam 
smearing is larger for the more distant Haro 11 we determine the 
velocity gradient is ~ 2.2 times steeper in Haro 11. The dynami¬ 
cal mass and density of a system scales as: 

M oc rv^ , p oc Al/r^ 


^ p oc (v/r)^ 

The morphological similarity suggests that inclinations are 
similar for Haro 11 and the Antennae. Therefore our results sug¬ 
gest that the mass density in the central parts of Haro 11 is 
2.2^ ~ 5 times higher than that of the Antennae. Of course, 
given the non equilibrium nature of mergers and the uncertainty 
of inclinations one should not take masses at face value, but the 
physical compactness of Haro 11 in any case argues for a higher 
density as compared to Antennae. 

The star formation rate density is known to , on average, scal e 


with the gas surface density as Esfr oc (jKennicutt 


19981. 


Since H = pr the difference in surface density is again 2.2 and 
hence Esfr should be a factor of 2.2^ '^ = 3 higher in Haro 11, if 
the gas mass fraction is similar in the two systems. 

The total star formation rate, SFR = r^SsFR, and we 
thus expect the total star formation rate of the Antennae to be 
2.2^/2.2^-^ = 1.6 times that of Haro 11. Hayes et a l. (2007) quote 
a probable SFR in the range 20 to 25 M© yr~^ and 


Madden et al. 


( |2014| ) finds 29 for Haro 11 whereas [Karl et ak] ( [2010 ) compile 


^ STScI News Release Number: STScI-2006-46, see: 
http://hubblesite.org/newscenter/archive/releases/2006/46 


various estimates for the Antennae in the range 10 to 20 M© yr“L 
The Ha luminosity (uncorrected for extinction) is about 3 times 
higher in Haro 11 and the IRAS 60yum luminosity is a factor of 
2 higher (based on the values in NED). Hence Haro 11 seems to 
be forming stars at 2-3 times the rate of the Antennae, contrary 
to the expectation from the mass scaling. 
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. Antennae (left, credit STScI and B. Whitmore) and Haro 11 (right, 


(20101, credit ES0/ESA/Hubble and NASA) 


put next to each other. The Antennae image has been rotated to the same apparent orientation as Haro 11 (which is shown as N up 
E left) and the images have been scaled so that the two systems look equally big on screen, for ease of comparison. The physical 
scale for the image of each system is shown by the red/green bar. In the image of Haro 11 we have labelled the three main knots 
according to their common designation A, B and C (Hayes et al. 2007). In the image of the Antenna we point out with arrows and 
roman numerals some regions of similarity which are discussed in text. In the lower right inset we show the image of Haro 11 at the 
same physical scale of the Antennae. The image of Haro 11 is a composite of observations from HST/ACS (filters E435W, E550W, 
E814W) an d VLT/NACO (K^). The inset plot in red (lower left) show the position-velocity (PV) diagram for the Antennae ( |Amram 


et al.|1992] ) and the green inset to the right shows the same for Haro 11 (this work). The PV diagrams have been extracted through 


the two nuclei in the Antennae and the corresponding for Haro 11 as seen in Pig.|^ i.e. knots C and Both PV diagrams cover the 
same total range in velocity (Av = 250 km s“^). The Antennae VP extend over 130"(16 kpc) while the one in Haro 11 covers 20"(8 
kpc) and the same velocity scale. As can be seen the shape and amplitude of the PV diagrams are remarkably similar. 


To get a proxy for the stellar mass ratio of the two galaxy 
systems we consider the K-band luminosities. NED gives Ks = 
7.2 for the Antennae, or (f or m - M = 31.9) Mk = - 24.7, while 
Haro 11 has Mk = -22.9 ( Bergvall & Qstlin||2002| . Hence the 
K-band luminosity of the Antennae is 5 times higher than that 
of Haro 11, and the stellar mass likely at least 10 times higher if 
one accounting for the fact that Antennae is more dusty and that 
Haro 11 due to its a stronger starburst should have lower MjLx. 
Hence, the specific star formation rate is at least 20 times higher 
in Haro 11 than in the Antennae. 

The formation of tidal tails in mergers have been investigated 
in several simulation studies. Early claims that lack of tails in¬ 
dicated very high dark matter fractions ([Dubinski et al.||1996| ) 
were countered by Springel & White) ( 1999) who showed that 


tails could be produced for high dark matter fractions provided 
the halo spin is sufficiently high. Hence the progenitor galaxies 
of the Haro 11 merger may have had smaller angular momen¬ 
tum or less extended dark haloes than the Antennae progeni¬ 
tors, both options suggesting more compact participating galax¬ 
ies in Haro 11. A retrograde merger would also lead to weaker 
tails, but this explanation is unlikely in the case of Haro 11 given 
the close similarity with the Antennae, which has been success¬ 
fully modelled with a pro-grade encounter (e.g. Karl et al.|2010| ). 
Certainly, dedicated hydrodynamical simulations would be key 
to better understand Haro 11. 


Other merging systems do not possess such obvious appar¬ 
ent similarity as Haro 11 and Antennae. Arp 299 has a somewhat 
similar morphology. It is halfway in size between Haro 11 and 
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Antennae and has a FIR luminosity 0.4 dex higher than Haro 11. 
Hence the SFR surface densities in Haro 11 and Arp 299 are 
comparable. While there is an abundance of gas in Haro 11 de- 
tect ed through various tracers it is remarkably devoid of cold 
Hi ( [Bergvall et al.||2000| |MacHattie et al.|Ml^. The SFR per 


unit mass is very high in Haro 11 and |Qstlin et al.| ( |1999| ) esti¬ 
mated that the mass of ionised gas in Haro 11 could be as high 
as lO^M©. The detection of [OIII] lines so far out suggest that 
the galaxy may be density bounded which would be consistent 
with the claimed direct detection of Lyman continuum radiation 


( [Bergvall et al. [200^ [Grimes et al.|2QQ7}|Leitet et al.|201l] ) 


To conclude, the morphological and kinematical properties 
of Haro 11 and the Antennae are stunningly similar, width the 
exception that the physical size of Haro 11 is a factor of ~ 4 
smaller and the implied mass density roughly 5 times higher. 
The star formation rate being twice as high in Haro 11 implies 
that the star formation efficiency is higher in Haro 11. 


5.6. On the origin of the line width 


It is known since long that the ionised gas follow a relation be- 
tween the line width cr and t he emission line strength of e.g. HJS 
( [Terlevich & Melnick|198l| ). The reason for such a relation has 
been discussed over the years, and while its discoverers have ar¬ 
gued that it originates in virial motions, this has been challe nged 
by others ( Green et al.|2010[[Moiseev & Lozinskaya|2012[ ) who 
argue that the relation with the emission line luminosity, hence 
star formation rate, indicates that star formation itself - through 
feedback - broadens the line. In the current context, this question 
is important as we have used the line width to estimate masses. 
Feedback is of course present also in Haro 11 and similar lumi¬ 
nous BCGs but the question is if feedback can produce lines with 
(T ~ 80 km s“^ i.e. affecting the bulk motions of the ionised gas 
rather than just producing broad low intensity wings? Any mat¬ 
ter in a galaxy must respond to the gravitational potential and 
in non rotating systems, the line width will always have a com¬ 
ponent from virial motions. In Haro 11 it is evident from this 
work, that the broad lines may be composed of many discrete 
velocity components and what could be wrongly interpreted as 
a very turbulent ISM is in fact a collection of moderately broad 
lines, which it takes high spectral and spatial resolution plus high 
signal to noise to reveal. The BCGs that have so far have their 
stellar motions investigated through Calcium triplet observations 
show that the Hii an d stellar line widths are remarkably simi¬ 
lar on a global scale (Ostlin et al. 2004 1 Marquart et al.[[2007 


[Cumming et al.|2008 this work). The motions of these predom¬ 

inantly young stars are not affected by feedback and they must 
have been born with this velocity dispersion. Hence, also the 
ionised gas velocity dispersion must be dominated by graviation. 

Studies of the evolution of the emission line width cr and 
Vrot!O' show both to increase with increasing redshift (e.g. 
Wisnioski et al. MTS] ). This can be intepreted as due to increased 


gas mass fractions, requireing higher cr to be stable according to 
the Toomre criterion ( [Toomre[|l964[ ). Therefore we must con¬ 
clude that the H ii line width seen in Haro 11 and similar systems 
to a large extent must be of gravitational origin. Numerical sim¬ 
ulations (e.g. [Bournaud et al.[[MT^ also show that turbulence 
drives star formation and not vice versa. The lines are not broad 
because of feedback from strong star formation, it is the turbu¬ 
lent medium which (by increasing the pressure and the Jeans 
mass) increases the star formation rate. 


6. Conclusions 

We have investigated the kimematics of ionised gas and stars in 
Haro 11, a luminous blue compact galaxy. Haro 11 hosts a very 
strong starburst in terms of specific star formation rate (or equiv¬ 
alently, birthrate 5-parameter), and is also one of the most nearby 
galaxies which matches z ~ 3 Lyman Break Galaxies in terms of 
far UV luminosity and surface brightness. The kinematical data 
have been collected with the ESO VLT and 3.6m telescopes. We 
also make use of imaging data from the Hubble Space Telescope. 

We find that the stellar and ionised gas kinematics to first 
order agree very well, both showing the same large scale fea¬ 
tures and similar velocity dispersions. Hence, wherever we can 
trace their motions, stars are moving with similar directions and 
amplitude as the ionised gas. This implies that the bulk of the 
ionised gas emission arises from gas whose velocities are gov¬ 
erned by gravitational motions. The irregularities in the velocity 
field must hence be attributed to real dynamical disturbances, 
providing further evidence that Haro 11 is a bona fide merger. 
This interpretation is also supported by the morphology, ex¬ 
emplified by the deep HST/NICMOS image showing telltale 
merger signs. The similar velocity dispersions in stars and gas 
also suggests that the Ha line width is primarily caused by virial 
motions. This suggests that emission lines provide a good probe 
of the kinematical state of dwarf starbursts, near and far. 

When analysing the ionised gas velocity field at higher spec¬ 
tral resolution with integral field spectroscopy we find several 
examples of double and triple line profiles. Only with high spec¬ 
tral resolution can these components be resolved. We find high 
and low velocity components that make however a small contri¬ 
bution to the total intensity and which would therefore be missed 
by low resolution spectroscopy. We used continuity arguments 
to associate the multiple line profiles with kinematical compo¬ 
nents. Notably, we find an elongated, kinematically cold, red- 
shifted arm. 

Using an off-centre IFU observation, we can trace the ionised 
gas ([Oiii] and Hy6) kinematics to large radii (r > lOkpc), where 
we find a fiat velocity curve offset from the centre with just 50 
kms“^. The line width is the same as in the centre. We discuss 
various alternative explanations for the outer velocity field, such 
as outflows and turbulence, but favour an interpretation where 
the outer halo is slowly rotating but mainly supported by velocity 
dispersion, implying a total dynamical mass of lO^^M©. 

We Anally compare Haro 11 to the famous Antennae system 
and And stunning similarities. Haro 11 appears to be a more com¬ 
pact version of the Antenna but otherwise showing very similar 
central morphology and kinematics. The Antennae is 10 times 
more massive than Haro 11 but only form stars at half the rate. 
One important difference is that Haro 11 lacks extended tidal 
arms, which may be due to more compact (denser) galaxies par¬ 
ticipating in the merger. 
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